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In Phys. Rev. A 71, 060312(R) (2005) the robustness of the local adiabatic quantum search to 
decoherence in the instantaneous eigenbasis of the search Hamiltonian was examined. We expand 
this analysis to include the case of the global adiabatic quantum search. As in the case of the 
local search the asymptotic time complexity for the global search is the same as for the ideal closed 
case, as long as the Hamiltonian dynamics is present. In the case of pure decoherence, where 
the environment monitors the search Hamiltonian, we find that the time complexity of the global 
quantum adiabatic search scales like iV 3 / 2 , where N is the list length. We moreover extend the 
analysis to include success probabilities p < 1 and prove bounds on the run time with the same 
scaling as in the conditions for the p — > 1 limit. We supplement the analytical results by numerical 
simulations of the global and local search. 
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I. INTRODUCTION 

Quantum adiabatic algorithms, first conceived in Ref. 
[2, have been developed as an alternative to the tradi- 
tional circuit model of quantum computation. An im- 
portant aspect of the adiabatic approach is its expected 
resilience to various kinds of open system effects, which 
makes it a promising candidate for robust quantum com- 
putation. The basic reason for this appealing feature is 
that the adiabatic quantum computer operates near the 
instantaneous ground state of a time-dependent Hamil- 
tonian and can therefore be expected to be insensitive 
to relaxation and open system effects among the excited 
states. The robustness of adiabatic quantum compu- 
tation has been discussed in the literature for different 
kinds of perturbations, such as relaxation Q, noise 0, 
decoherence 0, 0, HI , and unitary control errors . 

In this paper, we elaborate on the robustness of the 
quantum adiabatic search 0, Q, 13 m the presence of de- 
coherence in the instantaneous eigenbasis of the search 
Hamiltonian. In_particular,we extend the analysis of the 
local search 0, Q in Ref. y| to global adiabatic search 
. We further address the dependence of the asymptotic 
time complexity of the search procedure on the success 
probability for finding the searched item. 

To model adiabatic search in the presence of decoher- 
ence we consider master equations of the form (H = 1 
from now on) 

Q=-iA[H, Q ]-B[W,[W,Q]], (1) 

where g is the density operator, and A, B are constants 
such that B > 0. For convenience, we also assume 
that A > 0. The operator H is taken to be the time- 
dependent Hamiltonian of the adiabatic search 0, Q, 
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and W is a time-dependent Hermitian operator such that 
[H, W] — 0, which gives decoherence in the instantaneous 
eigenbasis of H. 

Decoherence with respect to the energy eigenbasis of 
a system has been shown to occur in certain regimes 
of system-environment interaction Apart from this 
general regime of weak interaction and dominating self- 
Hamiltonians one can find other settings where this type 
of decoherence occurs. One example is control errors in 
terms of fluctuations in the energy levels of the system 
[lfj . Note also that decoherence in the energy eigenbasis 
does not necessarily imply the particular form of Eq. Q . 
However, this choice provides a sufficiently simple model 
to be treated analytically, and yet results in non-trivial 
behavior. Moreover, the form of Eq. quite generally 
guarantees the solution to be a proper density operator 
at all times. 

Our results show that for A and list length N, it 
is sufficient with a run time of the order N in the global 
case and VN in the local case, as is the case for the ideal 
closed evolution. These results suggest that the adiabatic 
quantum computer is protected against this type of de- 
coherence, which is promising for phy sical realizations. 
In the wide-open case where A = jlJj , with the special 
choice W — H, the scaling changes to N 3 / 2 and N, for 
the global and local search, respectively. 

Although the concept of adiabaticity has a well de- 
fined meaning for closed systems, there is no obviously 
unique generalization to the case of open systems. One 
may consider several different generalizations focusing on 
various aspects of the standard adiabaticity. For exam- 
ple, in Refs. 0, 0, Q] the concept of adiabaticity is 
based upon Jordan block decompositions of the superop- 
erator that describe the evolution of the open system. A 
different approach, designed for systems that are weakly 
open, has been put forward in Ref. [l5j |. This latter ap- 
proach is based upon taking decoupling of the energy 
subspaces of the ideal Hamiltonian as the starting point. 
Here, we use a third method based upon a feature of the 
adiabatic quantum computer, namely that the evolution 
takes place near an eigenstate of the ideal Hamiltonian. 
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This property is in contrast with, e.g., implementations of 
holonomic quantum gates 0], where it is essential that 
the gate can operate on arbitrary superpositions, with- 
out too large errors. For the functioning of the adiabatic 
quantum computer, however, it is sufficient that the state 
of the system remains close to the correct eigenstate. In 
the present interpretation of the word "adiabaticity" we 
only require that the probabilities of finding the system in 
the instantaneous eigenstates of H (s) are conserved, but 
we allow superpositions of the instantaneous eigenstates 
to decay into mixtures. 

The structure of the article is as follows. In Sec. [H] we 
derive an integral equation from the original master equa- 
tion. Section IIIII gives a brief overview of the adiabatic 
search algorithm, and we calculate some quantities which 
are used later. In Sec. II VI we prove a sufficient condition 
for the success probability of the adiabatic search to ap- 
proach unity in the semi-open case. We consider both 
global and local search. Furthermore, a condition on the 
decoherence term of the master equation, which is used in 
the proofs in Sec. IIVI is examined. In Sec. we demon- 
strate that the wide-open case is essentially different from 
the semi-open case. We prove sufficient conditions for the 
success probability to approach unity in both the global, 
as well as the local case. We further prove that the suf- 
ficient conditions are also necessary. In Sec. IVII bounds 
on the run time for fixed success probabilities less that 
unity, are given. We end with the conclusions. 

II. THE MASTER EQUATION 

In this section we develop the basic equation used in 
our analysis. Since it is our intention to use adiabatic- 
ity we consider a fixed family of Hamiltonians H(s), 
parametrized, in a sufficiently smooth way, by the param- 
eter s. We assume that the parameter s is proportional 
to the time t, such that s — at time t — and s = 1 at 
time t = T. The time T is referred to as the "run time" . 
We also assume a family of Hermitian operators W(s), 
to be used in the master equation 

j t g(t)=-iA[H(t/T), Q (t)}-B[W{t/T),[W(t/T),g{t)}]. 

(2) 

We assume that both H(t/T) and W{t/T) are Hermi- 
tian and non-degenerate at each t. Moreover, we assume 
[H(t/T), W(t/T)} = for all t. By a change of variables 
s = t/T in Eq. J3J) one obtains 

-f p(s) = -iTA[H(s),p(s)] - TB[W{s), [W(s), p(s)}], 
as 

(3) 

where p(s) = g(sT). Let {\E n (s))} n be instantaneous 
orthonormal eigenvectors of H(s), i.e., H(s)\E n (s)) = 
E n (s)\E n (s)) . Since H(s) is assumed to be non- 
degenerate the arbitrariness in the choice of eigen- 
vectors is limited to phase factors. Let w n (s) de- 
note the eigenvalues of W(s). If we insert p(s) = 



J2 nm Pnm(s)\E n (s))(E m (s)\ into Eq. ©, the equations 
for the components p nm {s) become 

-rPnmis) = -iTA{E n (s) — E m (s)}p nm (s) 
as 

-TB{w n (s) - w m (s)} 2 p nm (s) 
-i[Z(s),p(s)] nm , (4) 

where Z(s) is an Hermitian matrix with elements 

^ urn 

(s)=i(E n (s)\E m (s)) (5) 

and p(s) is the matrix with elements p n m(s). 
Consider an ordinary differential equation 

— x(s) = Li(s)x(s) +L 2 (s)x(s), (6) 
as 

where Li(s) and L2{s) are s-dependent linear operators 
acting on the vector x. If the space is finite-dimensional, 
Li(s) and £2(5) are sufficiently well behaved as functions 
of s, and if £i(s) fulfills [Li(s), Li(s')] = for all s, s' , 
then Eq. © can be rewritten as the following integral 
equation 

x(s) = efo L i( s ') ds 'x(0) 

+e I " Li( s ')ds' r e -/o 8 'ii(«")^" j L 2 ( s ')- E ( s ') ds '. 

(7) 

In the present case we have x(s) — p(s) and 

{L 1 (s)p{s)) nm = -iTA{E n (s) - E m (s)}p nm (s) 

-TB{w n (s) - w m (s)} 2 p nm (s), (8) 
L 2 (s)p(s) = -i[Z(s),p(s)}. (9) 

Hence, Eq. J2J| can be rewritten as 

Pnm(^) — ^ ^ ^ Q ( ^Pnm(O) 

_ ie -iTAR nm {s)-TBQ nm (s) 

x /' S e ™„ m ( s ') + T S Q„ m (,') [Z(s0ip(s0]nmds ' i 
JO 

(10) 

where 

Qnm(s)= [ [w n (s')~W m ( S ')} 2 ds', (11) 
JO 

Rnm{s)= I [E n { S ') - E m { S ')]ds' . (12) 

Jo 

III. THE SEARCH PROBLEM 

Our task is to find a single marked item in a collection 
of N items. The iV-element search problem is associ- 
ated with an iV-dimensional Hilbert space with orthonor- 
mal basis {\k)}k=i, where the marked item corresponds 



3 



to e {1^)1^1^ We consider the following family of 
Hamiltonians PUT], 



H(t/T) = -(1 - t/T)\i/,)W - t/T\»){n\, 



where 



1 N 



(13) 



(14) 



The ground state of the initial Hamiltonian H(0) is 
easy to prepare, while the ground state of H(l) gives the 
solution of the search problem. Hence, if we assume the 
evolution is adiabatic, this family of Hamiltonians takes 
us from the initial state \ip) to the marked state \ fi), and 
thus solves the search problem. A crucial observation is 
that the problem is essentially two-dimensional, since the 
space Hr spanned by \ip) and is the only relevant sub- 
space. We may Gram-Schmidt orthogonalize {\ip), \n)}, 
yielding the orthonormal basis 



N\H) ~ W 



(15) 



We represent the Hamiltonian family in Eq. (|13|) . re- 
stricted to the space TLr, as matrices in the basis defined 
in Eq. |TH)l. After the substitution s = t/T, we get 



M(s) 



,N-1 




N 
, N-1 
'"AT 



which has the eigenvalues 



and orthonormal eigenvectors 
,(0) _ \/2 



(16) 



(17) 



where 



VA 2 - AD 
VA 2 + AD 



A(s) = E 1 (s)-E (s) 



L N 




4 0) (,) 


a(A + Dy 

L jv - 




4 1} ( S ). 



(18) 



l + (iV-l)(2s-l) 2 



D(s) 



-1 + 2s 



N 
N - 1 
N 



(19) 
(20) 



The two relevant eigenvectors of H (s) are 



\Eo(s)) 



\Ei(s)) = e?\sM 



o(°) 
-i 

(i)/ 



WW 



(21) 



where \Eq(s)) is the ground state. 

We next calculate the matrix Z(s), defined in Eq. JSJ. 
In order to save some efforts we first notice that Zqo(s) 
as well as Zn(s) are identically zero, since eg°\ e£ , e^ 1 , 

and are real valued. Secondly Z(s) has to be an 
Hermitian matrix, which follows from differentiation of 
the relation (E n (s)\E m (s)) = 8 nm . 

By combining Eq. JSJ) with Eqs. ijTSj) . ijlilj) . and J3UJ, 

one obtains the following result 



Zm (s) = 



-WN - 1 



1 + (N — 1) (2s - 1)' 



(22) 



Define 



Z(s) = \Z 01 (a)\ 



1 + (N- l)(2s - If 



(23) 



Note that the maximum y/N — 1 of Z(s) is obtained at 
s = 1/2. A useful property of Z is 



J*Z(s)ds<^, 



(24) 



for all N > 1. 

Concerning the reduction to the two-dimensional prob- 
lem, one should note that, due to the assumption 
[H(s),W(s)] = 0, the operator W(s) block-diagonalizes 
with respect to the eigenspaces of H(s). This implies 
that the solution of Eq. J5J) remains in the subspace TCr 
once started there. Thus, we only need to consider the 
master equation in Eq. (0 on this particular subspace. 

Define 



Q(s) = Qio(s), R(s)=R 10 (s), (25) 
Y(s) = poo(s) - pu(s) = 2p o(s) - 1. (26) 

From Eq. I|10(l one can derive the following equation for 
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Y{s) = F(0)-4 



-TBQ(s') y/i 



Z(s)Rc 



-iTAR(s') 



Pio(O) 



ds 



-4 / e - TB ^') / gT_BQ(s")j^ e 



e -iTAR(s') e iTAR(s") 



Z{s')Z{s")Y{s")ds"ds'. 



(27) 



As seen from Eq. (|26[) . the value of F(s) determines to 
what degree the system is in the ground state, i.e., Y = 1 
when the system is in the ground state, and Y = — 1 when 
the system is in the first excited state. 

IV. SEMI-OPEN SYSTEM 
A. Global search 



where if is a constant independent of N. This assump- 
tion essentially means that the fluctuations of T in s are 
not allowed to grow with the list-length N. The restric- 
tion to an TV-independent K can be relaxed. We discuss 
this issue, as well as the condition in general, in Sec. lIV C\ 
From Eq. it follows that 

|r(0)-y( S )|<4|J( S )|+4|7( S )|, (34) 

where 



We begin with a comment on our terminology re- 
garding the "openness" of quantum systems. With a 
"wide-open" [llj system we intend a system governed 
by Eq. Q), with A = 0. As a contrast to the wide-open 
case we define "semi-open" as A > 0. Finally, a "closed" 
system denotes the case B = 0. 

We define 



r(s) = Wi(s) - w Q (s), 

R(s) = f A(s')ds', 
Jo 

Q( S ) = f r 2 ( s ')<fe', 



5 = min A(s) = , , 

«e[o,i] V^V 

7 = min F(s). 

*e[o,x] 



(28) 
(29) 

(30) 

(31) 
(32) 



In the proofs we use the following restriction on the de- 
coherence term 



Z(s) 



ds 



r 2 (s) 



ds < K, 



(33) 



(35) 



I±(s) 
u±(s') 



e- T ^^ ±iAR ^Z(s')u ± (s')ds',(36) 
e T ^ BQ ^ ±iAR ^Z(s")Y(s")ds", 

(37) 



and 
J(s) 



-TBQ{s') 7 fJ 



Z(s')Re 



-iTAR(s') 



Pio(0) 



ds'. 



(38) 

For the sake of simplicity we separate the treatment of 
the terms I and J. 



1. The term I 

To begin with we define the following integral and per- 
form a partial integration 
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G±[g](s) = f e- T ^ s '^ AR ^Z( S ')g( S ')ds' 
Jo 



-T[BQ(s')±iAR(s')] 



ds' 



Z(s')g(s') 



-ds' 



. _ p -T[BQ(s)±iAR(s)] _ 

Z(s)g(s) 



i Vn^i 



9(s) 



- f p -T[BQ(s')±iAR( s ')]_^_ 
T 



s=0 



Z(s') 



ds' \ B %l(s')±iAj§( S ')^ 



g(s')ds' 



+ 1 [' c -T\BQ(s')±iAR(s')] Z ( S ') 



T 

If we put g(s) = u±(s) into Eq. the result is 

Z(s) 



B d S{ S ')±^(s')ds 



7 g(s')ds'. 



I ± (s) = G ± [u±](s 
1 



T[BQ(s)±iAR(s)] _ 



T 
1 

i 



e T[BQ( s ')±Mfl( s ')] Z(s / )y(s / )ds / 



-T[SQ(s')±iAfl.(s')]_ 



Z(s') 



(39) 



t7 B ^ (s0 ±a^( S ') 



ds' \ B ^(s')± l A^(s')J Jo 
Z(s')Y{s')ds'. 



,T[BQ{s")±iAR{s")} z ^ l)Y ^, )ds n ds , 



(40) 



r 



In the next step we calculate an upper bound to |i±(s)|. 
To do this we use that |^(s)| < 1. Moreover, since 
s' > s" it follows that Q(s') > Q(s"), which im- 
plies exp{-TB[Q(s') - Q(s")}} < 1. We also use that 
dQ/ds = T 2 (s), dR/ds = A(s), and Eq. to obtain 



\I±(s)\ < 



1 7T 

T 2 
1 



Z(s') 



ds' \BT 2 {s')±iAA(s') 
Z 2 (s') 



ds' 



T J Q x /B 2 T i (s')+A 2 A 2 {s') 

1 7T Z{s) 

T2 yfWfH^s) +^4 2 A 2 (s)' 



ds' 



Now we use the fact that Z{s) < y/N — 1, as well as 
Eqs. (|3*TJl and lj3~2l . which result in 



ds' \BT 2 {s')±iAA(s') 



To treat the integrand in Eq. 1|42J) . we define 



X( S ) = 



Z(a) 



ds \BT 2 (s)±iAA{s) 



which can be estimated as 
X(s) < ' 



B 



ds 



Z(s) 



+ 



B 2 1 i + A 2 S 2 

A 

' B 2 1 i + A 2 5 2 



Z(s) 
Z(s) 



r 2 (s) 

A(s 



(43) 



The function Z{s) is symmetric around s = 1/2 and in- 
creasing on the interval [0, 1/2]. It follows that the func- 
(41) tion \dZ/ds\ is also symmetric around s = 1/2, and that 
\dZ/ds\ = dZ/ds on [0, 1/2]. Thus 



ds' < 



d , , s 



ds' < 2\/N - 1. 



(44) 

dg i Since both Z{s) and A(s) are symmetric around s — 1/2, 

Z(s)\dA/ds\ has the same symmetry. Moreover, A(s) is 
increasing on the interval [1/2, 1]. Hence, Z(s)\dA/ds\ = 
(42) Z(s)dA/ds on [1/2, 1], which leads to 



Z(s) 



d a / \ 
^-A(s) 
ds 



ds < 



Z(s) 



d a / \ 
ds 



ds 



N-l „ 
< 2j < 2. 

TV ~ 



(45) 



6 



By combining Eq. Jj3J with Eqs. (gSJ, |g2J|, and l|15^. 

one obtains 



x W < 2V7 ^ ' BK + 2A 



(46) 



By combining Eqs. l(3T)l. g2J|, and lO it follows 

that 



1^01 < 77 7F 



BK + 2A 



2 T y ^B 2 ~f 4 N + A 2 B 2 ^N + A 2 J ' 

2. The term J 
Equation (|38[1 can be rewritten as 

J{s) = Re ( Pl0 (0) f e- T ^ s '^ AR ( s '»Z( S ')d S ' 



(47) 



(48) 

We write \J(s)\ < \pw(0)\ \G+[l](s)\, with G as in 
Eq. I|39|) with g(s) = 1. By a line of reasoning very 
similar to the one in the previous section one obtains 



|G+[1](*)| < 



1 

TJ 



Z(s') 



1 



ds' Vsr 2 (s') + iAA{s') 
1 



+ 



Ty/N y/B^iO) + A 2 

1 Z(s) 

T y /B 2 T A {s)+A 2 A 2 {s)' 



ds' 



(49) 



which results in 



\J(s)\ < |/oio(0)| 



1 



1 



Ty/N y/B 2 r 4 (0)+A 2 

pM0)\ N 3 



Ipio(0)| 



T y/B 2 ^N + A 2 
,N BK + 2A 
T B 2 ^N + A 2 ' 



(50) 



3. Collecting the results 
By combining Eqs. (|.j4p. (|47|) . and (JSUJl, one obtains 

|y(o)-y(a)| < 4|p! (o)| 1 



T^A 7 y/B 2 T\Q) + A 2 
4^10(0)1-^ " 



T y/B^N + A 2 
N BK + 2A 
T B 2 ^N + A 2 
N 1 



+4|pio(0)| 
+8tt 

+2tt 



T y/B 2 -fN + A 2 
N BK + 2A 
f B 2 -f 4 N + A 2 ' 



(51) 



Note that 

B 2 j 4 N + A 2 > A 2 , 
B 2 T 4 (0) + A 2 > A 2 , 

which, as long as A > 0, implies that 
Ipdo(0) - /»oo(s)| < 2jp X0 (0)] 



(52) 
(53) 



1 1 



tVn a 

t i N BK + 5 A 
+2|pio(0)|- 

N BK + 6A 



T A 2 



(54) 



where we have used Eq. (j26(l . The leading term in 
Eq. is proportional to AT/T. In the N/T — > oo limit 
the probability to find the system in the groundstate is 
conserved. Thus, it is a sufficient condition for adiabatic- 
ity that T > N. 



B. Local search 

As is well known [l(, the standard global adiabatic 
search algorithm performs as the classical search, but 
the local adiabatic search 0,0 outperforms the classical 
algorithms. Here we further elaborate and extend the 
analysis of Ref. [j| concerning the effect of decoherence 
on the local adiabatic search. The idea behind the local 
search is to adjust the speed of parameter change in such 
a way that more time is spent near the minimum energy 
gap. This is obtained by a reparametrization of the func- 
tions H{s) and W(s) to H(f(s)) and W(f(s)), where / 
is a strictly increasing, sufficiently smooth function, such 
that /(0) = and /(l) = 1. The new master equation 
becomes 



-p(s) = -iAT[H(f(s)),p(s)} 
as 

-BT[W(f(s)),[W(f(s)),p(s)}]. (55) 

By a change of variables r = /(s) and by defining p(r) = 
p(/ _1 (r)), Eq. ||SSJ becomes 

-fp(r) = -iAT^(r)[H(r),p(r)] 
dr dr 

-BT^(r)[W(r),[W(r),p(r)}}. (56) 

Note that p(0) = p(0) and p(l) = p(l). 

The equation for Y(r) — 2p ao (r) — 1 is fully analogous 
to Eq. (|27|) with Q and R replaced by Q and R, where 



Q(r) 



R(r) 



T 2 {r')^—dr', 
dr 

A(r')V—dr'. 
dr 



(57) 
(58) 
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The choice 



By using Eq. it follows that 



L 



1 



LJo A 2 (r') 
1 



dr', 



o AV) 



dr' 



(59) 
(60) 



gives_the optimal efficiency for the quantum search prob 



is the opt ii 
lem 0,0,01 ■ One obtains 



Q{r) = T 



LJo A 2 (r') 



dr' 



R(r) = \ 



1 



A(r') 



dr', 



£ = 



arctan(V^V - 1) < - 



N 



2 % /lV^T 



Define 



. r 2 ( s ) 

mm 

s e[o,i] A(s) 



(61) 
(62) 

(63) 
(64) 



As in the global case we let \Y(0) - Y(r)\ < 4|J(r)| + 
4|/(r)|, and as before we separate the treatment of the 
terms I and J. The terms I(r) and J(r) are defined as 
in Eqs. lpjS )l -l|3*7" |l and gSJ, but with Q and R replaced 
by Q and R, respectively. 



1. The term I 

The analogue of Eq. I|40|). is obtained by replacing 
Q and R by Q and R, respectively. We use dQ/dr — 
T 2 (r)/[LA 2 (r)], dR/dr = l/[LA(r)], exp{-TB[Q(r) - 
<2( r ')]} < 1 if r > r', as well as Eq. (J23J, which result in 



|/±MI < 



1 7T 

T2 
1 

+ T 



d 

d? 



Z(r>) 



L A 2 (r' 

Z 2 (r') 



L L A(r') 



d/ - ' 







1 7T 

+ T 2 

Xtt r 



L 

+ T 



R 2 1 r4(rQ 
-° L 7 A 4 (r') 

Z(r) 



=dr' 



4 2 1 . 1 



-° L 2 A 4 (r) 



A 2 - 



d 

dr 7 



L 2 A 2 (r) 

Z(r')A(r') 



A(r') 



dr' 



"T 2 



Z 2 (r')A(r') 

^a^ + ^ 2 
Z(r)A(r) 



=dr' 



D2 r*(r) 
■° A^(7) 



A 2 



(65) 



|/±(r)| < 



L 7T 

T 2 
L 



Z(r')A(r') 
dr' \ B M±a 



dr' 



Lit Z(r)A(r) 
T2 ^B\ + W 



Z 2 (r')A(r')dr' 



(66) 



One may use Z(r)A(r) < ^{N - l)/N and Eq. to 
obtain 



\I±(r)\ < 







-f 


T2 


f 

la 


dr' ^ 



Z(r')A(r') 



ttL IN -I 



A(r') 
1 



dr' 



# \/B 2 (^ 2 + A 2 



(67) 



The derivative contained in the integrand on the right- 
hand side of the above equation is performed by first 
multiplying both the numerator and the denominator by 
A(r') and then using Z(r)A 2 (r) = - l/N. More- 
over, we use Eq. I|64|l . The result is 



Z(r')A(r') 



< 



B 



B 2 ( 2 + A 2 J 

a r 1 



B 2 C+A 2 



Z(r') 
Z(r') 



dr' 



7 rV; 



dr' 
dr'. (68) 



Note that in the inequality above we have extended the 
integration interval from [0,r] to [0,1]. By combining 
Eqs. l|33 Jl . P ty . l(67 | . and ip | one obtains 



\I{r)\ < 



7T 2 y/N 



1 



+ 



2 T y^/B 2 C, 2 +A 2 B 2 ( 2 + A 2 
tt 2 VN 



N 



BK 



4 T V N - 1 B 2 C 2 + A 2 ' 



(69) 



2. The term J 



As in the case of global search one obtains |J(r)| < 
|pio(0)| |G+[l](r)|, with G as in Eq. with g(r) = 1, 
and Q and R replaced by Q and R. Note that p(0) = 
p(/ _1 (0)) = p(0) and thus pio(0) = pio(0). Using similar 
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arguments as in the previous sections, one obtains 

Z(r')A(r'[ 



|G+[l](r)| < - 





-f 


f 

10 


dr' \ 



A(r') 



dr' 



1 



+ 



T N yjB 2 T 4 (0) + A 2 
L Z(r)A{r) 



(70) 



Next we use Z(r )A(r) < y/(N - l)/N together with 
Eqs. tSJ|, and iJBBJl, yielding 



+ ^ lo(0 4 ^ 2 r4 + ^ 



7T ..VJV BAT / AT 
-2lPM0)| — B 2 C 2 + A 2Va73i 

7T V/V 2,4 

-2^10(0)1— B 2 C2+A 2- 



5. Collecting the results 



(71) 



We combine Eqs. (jSHJ and ((HI together with |Y"(0) - 
?(r)| < 4|7(r)| + 4| J{r)\, A > 0, B 2 C 2 + A 2 > A 2 , and 
B 2 r 4 (0) + A 2 > A 2 . Furthermore, we assume N > 2, 
which gives yjN/{N - 1) < a/2, and results in 



|poo(0)-A)o(r)| < (|pio(0)| 



-(3|pio(0)| + 2tt 

i~ , 7T 1 

-IpioWI^-. 



21 I 5 jT 

7T aZ/V 



4 T 



(72) 



The leading term of this equation is proportional to 
s/~N/T. Hence, a sufficient condition for adiabaticity is 
T > VN. Equation O can be expressed in terms of 
the solution p(s) of Eq. lfS1))l by the substitution r = f(s). 
However, this change of variables does not affect the con- 
dition for adiabaticity. 



C. The condition 

Both for the global and local search we have restricted 
the decoherence term of the master equation to fulfill the 
condition in Eq. (|33[) . This condition puts a limit on the 
fluctuations of T with respect to s, as a function of the 
list-length N. However, one may note that the condition 
can be relaxed, in the sense that we may allow K to be 
dependent on N. By inserting K = cN a , where c > 
and a > are constants, into Eqs. 1|54[) and (|72[1 . one 



finds that the new scaling becomes N 1+a and N x / 2+a 
for the global and local search, respectively. Thus, the 
scalings get worse, but one finds that there is a "window" 
< a < 1/2 where the local search still outperforms the 
classical search, although not with the optimal efficiency 
of the Grover search [l7| ■ 

In the following we attempt to obtain an indication of 
to what extent the condition in Eq. (|33fl is restrictive or 
not. For the sake of simplicity we do so for the more re- 
strictive case where K is independent of N. We assume 
r(s) = r/(A(s)), where r\ : (0,oo) — > (0, oo), is an in- 
creasing and sufficiently smooth function. Now, a larger 
value of i?r 2 (s) implies stronger decoherence. Thus, we 
have chosen to consider a decoherence model for which 
the strength of the decoherence increases with the energy 
difference between the two instantaneous eigenstates, but 
we do not specify any details on how it grows. We intend 
to derive a sufficient condition on rj for T to fulfill the 
condition in Eq. (|33|1 . 

By the assumptions on rj and the form of A, it follows 
that r is decreasing on the interval [0, 1/2] and increasing 
on [1/2, 1]. Partial integration implies 



Z(s) 



d r, , , 



r 1 d 

ds = 2 Z(s)—ri 2 (A(s))ds 

J 1/2 ds 



-2^/lsH^Tr, 2 (N- 1/2 ) 



r, 2 (A(s)) — Z(s)ds 



1/2 



ds 



(73) 



We use rfiN- 1 / 2 ) > and 2rf{\ yN - l/N < 2r) 2 (l). 
Moreover, we use Z(s)A 2 (s) = \/N — l/N, which leads 
to 



Z(s) 



ds 



T 2 (s) 



N 



1/2 



A 3 (s) ds 



(74) 



Only the integral on the right-hand side of the above 
expression may potentially grow when N increases. If we 
treat this term separately we find 



N J 1/2 A 3 (s) ds 3 ~ sfN 



l/VN 



r^A) 
A 3 



dA. 



This implies that a sufficient condition for the inequality 
in Eq. (|33fl to hold, is that the expression on the right- 
hand side of the above equation, is bounded. From the 
assumptions on rj it follows that the limiting behavior of 
this integral is determined by the asymptotic behavior of 
the function rj as A — > 0. Assume that r](x) < Cx a for 
some constants C and a > 0. Then it follows that 



1 r, 2 {x) 



1 



Vn J-L. x 3 



dx < 



C 2 InN 
2 y/N 

C 
2(a- 



AT2- 



a = l, 
(75) 
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In the case where a = 1, the function T satisfies the in- 
equality in Eq. ||33J, since (\nN)/\/N is a bounded func- 
tion on the nonzero natural numbers. In the case where 
(7^1, one can see that the right-hand side of Eq. (|T5|) 
contains two terms. None of these terms increase in mag- 
nitude, with increasing N, if a > 1/2. Hence, we may 
conclude that for all a > 1/2 the condition in Eq. I(33|) is 
fulfilled. In other words, for a sufficiently smooth, pos- 
itive, and increasing function 77, the only condition we 
have required is rj(x) < Cx a with a > 1/2. 



V. THE WIDE-OPEN CASE WITH T(s) = A CT (s) 

The wide-open case is obtained if we let A = (and 
B = 1 for convenience) which in the global search case 
yields the master equation 



it follows that Q(s') - Q(s") > N~ a (s' - s"). This can 
be used, together with \Y(s)\ < 1, to show that Eq. l(77|) 
leads to 

\Y(p)-Y(s)\ < 4 f /" ' e -TN-'(s'-s") Z (s')Z{s")ds"ds' 
Jo Jo 

+4jpio(0)| C e' TN ~° s ' Z{s')d S '. (79) 



One may use Z(s) < y/N — 1, integrate, and use Eq. flty 
to obtain 



|poo(0) - poo(*)| <(2|pio(0)|+7r): 



T 



(80) 



Thus, we have proved that a sufficient condition for poo(s) 
to approach Poo (0) is T > N a+1 / 2 , if a > 0. 



-p(s) = -T[W(s),[W(s),p(s)}] (76) 
and in the local search case yields 



[W(r),[W(r),p(r)]], (77) 



where we have changed variables to r = f(s) as described 
in Sec. ITVRl 

Even with the assumption that W(s) should be simul- 
taneously diagonalizable with H(s), we still have a rather 
large freedom in the choice of W, in terms of how T de- 
pends on the parameter s and the list-length N. In the 
previous sections it was only necessary to require that 
r should fulfill the condition in Eq. (|33l) . In this case, 
however, we need to specify more about W. 

Our main object of study in this section is the case 
W(s) — H(s), which corresponds to T(s) — A(s). How- 
ever, in order to prove that in the wide-open case the 
asymptotic behavior is dependent on the choice of W, 
we consider the generalization T(s) = A CT (s). Note 
in particular that if a = n, where n is an odd natu- 
ral number, then T(s) = A™(s) can be obtained with 
W(s) = 2 n - 1 (H(s) + (1/2)1)™. 

Here we prove that in the wide-open case, with T(s) = 
A CT (s), the global and local conditions for adiabaticity are 
T > N a+1 / 2 and T > N a , respectively. We first prove 
that these are sufficient conditions. Thereafter follows a 
proof of necessity. 



A. Global search: Sufficiency 

Let A = and B = 1 in Eq. I|27(l and assume that 
r(s) = A CT (s), which gives 



Q(s)= /" A 2 %s')d.s'. 
Jo 



(78) 



Assume a > 0, and use A(s) > S = 1/y/N. It follows 
that Q(s) > N' a s, for all s > 0. Moreover, since s' > s", 



B. Local search: Sufficiency 

Here we require a to fulfill a > 1. As before we let 
A = and B = 1 in Eq. J57J, but we replace Q with Q 
defined as 





1 r 



LJo 



A 2a - 2 (r')dr' 



> 



LN° 



(81) 



with L as in Eq. 1)63(1 . and where we in the inequality 
have used A(s) > 8 = 1/vN and the assumption a > 1. 
Inserting Eq. (jSTJ into Eq. <j77|l yields 



\Y(0)-Y(r)\ 



< 4 



^0 



e-^'- r "^ LN ^Z(r')Z(r")dr"dr' 



+4|p 10 (0)| / e-^ LN ^Z(r')dr'. 



(82) 



By a reasoning very similar to the derivation leading from 
Eq. J7HJ} to Eq. (gQJ it follows that 



|F(0)-Y"(r)| < 2tt 



LN a - l x/W~^l 
T 



^IP.oWI ^-'r 31 . (83) 



T 



which, together with Eq. (|63() . gives 



7T AT 17 

|poo(0) - poo(r)\ < (2|pio(0)| + tt)- — . (84) 
Thus, a sufficient condition for adiabaticity is T ^> N 17 . 
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C. Necessity 

In Sees. IV Al and lV Bl we proved sufficient conditions for 
conservation of the probability to find the system in the 
ground state, throughout the evolution. Here we prove 
that these sufficient conditions are also necessary, but 
under the simplifying assumption that we start in the 
ground state of the system. Moreover, we focus on the 
end point s = 1, i.e., we consider the success probability 
V = Poo(l), and necessary conditions for p to approach 
unity. 

Due to the effectively two-dimensional Hilbert space of 
the problem we may use a Bloch vector representation in 
terms of which the state of the system and the Hermitian 
operator W read 



p(s) = -[l +v (s)-a], 



(85) 



Wis) 



K[tflo(«)+wi(«)]i + u(s).<r}, (86) 



where v and u are elements in R 3 , ||u|| < 1, and er = 
{o~ x , o~ y , o~ z ) are the Pauli spin operators. It follows that 
u(s) — ~T(s)q(s), where q(s) is the Bloch unit-vector 
representing the ground state of H(s). 

If we use the assumption T(s) = A CT (s), then Eq. i|76|) 
can be rewritten in the Bloch representation as 



where 



v = -TZ(s){v-q{s)[q(s)-v]}, 



£(s) = A 2 *(s), *>0. 



(87) 



In the local case Eq. I|77|) again becomes Eq. 18 7|) . but 
with the parameter s replaced by r and 



£(r) = A 2 » 



dr 



a > 1. 



(89) 



Most of the material presented in this subsection is in- 
dependent of whether we consider the global or the local 
search. The difference is exclusively carried by the func- 
tion £. Due to this we will in the rest of this subsection 
use the variable s, although we should, to be correct, use 
r in the local search case. In Sec. IV (j"2"l we properly use 
the variable r again. 

Let and \tp) defined in Eq. I|15fl correspond to the 
north and south pole of the Bloch sphere, respectively. 
The Bloch vector q, corresponding to the instantaneous 
ground state, can then be expressed as 



i 



A(s) 



y/N - 1 N - 1 

2s^—rz ,0,1 -2s 



TV 



N 



(90) 



Note that q moves in the x-z plane only, and that q(Q) = 
(0,0,1). Let 8(s) be the angle between q(s) and q(0). 
This angle can be expressed as 



cos(0(s)) = q(s) ■ g(0) 



1 - (TV — l)(2s — 1) 

]V v /l + (/V-l)(2s-l) 2 ' 
(91) 



For N > 2 the right-hand side of the above equation is 
strictly decreasing on the interval < s < 1, beginning 
at the value 1 and ending at the value 2/N — 1. It follows 
that we can choose 8(s) as an increasing function with 
range in [0, n]. 

Since the system is started in the ground state p(0) = 
\E o (0))(E o (0)\, the initial condition is v (0) = q(0). Since 
both q(s) and v(0) are fixed to the x-z plane it follows, 
from the form of Eq. I|87() , that the motion of v is fixed 
to the x-z plane. Thus, we may restrict our analysis to 
this plane only. 

Lemma 1 Let qo, q' , and q" be three unit vectors in IR 2 , 
and let v be such that \\v\\ < 1. Suppose that the angle 
between qo and q" is greater than or equal to the angle 
between qo and q' . Moreover, let both these angles be in 
the interval [0, tt] . Then 



implies 



v q >qo <l 



q">qoq"- 



(92) 



(93) 



Proof. Consider a coordinate system such that qo 
corresponds to (1,0), q' to (cos#i,sin#i), and q" to 
(cos 02, sin 02 ), with < 6\ < 62 < 7r. In the cases when 
61,62 assumes possible combinations of 0,7r, the lemma 
is trivially true. Hence, we assume < 6\ < 62 < tt. 

If we let v — {x,y), then Eq. becomes cccos#i + 
ysm6\ > cos 8%. By using sin#i > 0, this can be rewrit- 
ten as 



y > (1 — x) cot 81 



(94) 



Since cotangent is a decreasing function on the interval 
(0, 7r), and since 1 — x > 0, it follows that the right-hand 
side of Eq. I|94|) is larger than or equal to (1 — x) cot 62- 
The above procedure can be reversed to obtain Eq. 19311 . 
□ 



Lemma 2 Ifv(s) is the solution of Eq. j6'7| ) with initial 
condition v(0) = q(0), and with q(s) as in Eq. 1^90)) . then 



v(s) ■ q(s) > q(0) ■ q(s), VsG[0,l]. 



(95) 



Proof. The strategy of this proof is to make a Cauchy- 
Euler polygon approximation of Eq. (|87|) . and show an 
inequality similar to Eq. I|95|l for every step size. Since 
the solution v(s) is obtained in the limit of small step 
sizes, the solution will satisfy the inequality in Eq. (|95ll . 

Let (sk)kLo be a partition of the interval [0, 1] such 
that so = 0. The partition has the step size Ds = 1/M. 
Let qk = q(sk) and let 



v 

Vk+l 



qo. 



v k - DsT£(s k )[v k - q k (q k ■ v k )]. (96) 
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Let v^^s) be the linear interpolation of Vq, v\, . . . , Vk, 
meaning that on the interval [sfe,Sfe + i] the function 
v^ Ds ^(s) is defined by 



(s) = v k+1 



Sk+l - s k 



Vk- 



Sk+1 



Sk+l 



Sk 



(97) 



Now we prove that the approximation v^^s) con- 
verges uniformly on the interval [0, 1] to the solution v(s), 
as Ds —> 0. Define 



f( S) v)=-n(s){v-q(s)[q(s)-v}}. 



(98) 



In both the global and the local case the function / can 
be shown to be continuous on I x I 3 , and to have a 
continuous partial derivative with respect to s. If one 
uses the assumption a > in the global case, and a > 1 
as well as the additional assumption N > 2 in the local 
case, one finds 

< £(s) < 1. (99) 

It follows that / fulfills the Lipschitz condition 

\\f(s,x)-f(s,y)\\<T\\x-y\\, (s,x),(s,y) e [0,l]xM 3 

(100) 

If we assume that Ds < 1/T, then Eq. I|99|l implies 

\X-DsT$(s)\ < 1, Vse[0,l]. (101) 

By induction using Eqs. i|96|) and l|101|) one can show that 
||ffc|| < 1- This in turn implies that the linear interpola- 
tion v^(s) fulfills 



\v (Ds) (s)\\ < 1, V*€[0,1]. 



(102) 



The inequality \\vk\\ < 1, together with Eq. I|9t)|) . can be 
used to show that 

||wfc+i - wjtll < DsT^a k )\\v k \\ < (s k +i - s k )T. (103) 
This gives 

\\vW{8")-v^^)\\<T\^-8'\, s',s"e[0,l]. 

(104) 

The fact that the partial derivative of /, with respect 
to s, is continuous on 1 x I 3 , implies that there exists a 
constant C such that 

\\f(s",v)-f(s',v)\\ <C\s"-s'\, (105) 

for all s',s" <E [0,1], and all v € R 3 : ||«|| < 1. By 
combining Eqs. i|100|) - (|lU5|) one obtains 

\\f(s",v^\s")) - f(s',v^\s'))\\ < (T 2 + C)\s" - s'\, 

(106) 

for all s' , s" € \0, !]■ This can be used to show that 



^(s")-^)^')- / /( S ,« (a) ( S ))d s 



(107) 



for all s', s" S [0, 1], and where e = Ds(T 2 + C). Since 
both / and ' are continuous, Eq. 11071) implies that 
v^ ^ is an e-approximative solution [19ll20| to Eq. JS7J). 

Since / is continuous and fulfills the Lipschitz condi- 
tion in Eq. (|lT)Ujl it follows [H HJ from Eq. ifTUTjl that 



\v W {s)-v(s)\\ < \\v (Ds \0)~ v(0)\\e 



Ts 



-Ds 



(T 2 + C),ts 



T 



(108) 



for all s G [0,1]. Since v(0) = »< ft >(0) = q(0) it fol- 
lows from Eq. (fTng|) that v^(s) converges uniformly on 
[0, 1] to the solution v(s), as Ds — > 0. (One may note 
that the solution v(s) exists and is unique, since / is 



continuous and Lipschitz on [0, 1] 



.) The solution 



v(s) is continuous as it is the limit function of the con- 
tinuous and uniformly converging functions v^ a "Hs) on 
the subset [0, 1] of K. (See Theorem 7.12 in Ref. f2lj.> 

In the following we focus on the proof of the inequality 
in Eq. I|95|l . We use the fact that the vectors v k and qk 
are restricted to the x-z plane of the Bloch sphere. By 
induction we prove that Vk ■ qk > <7o ' 9fe- 

By using Eq. I|96|l one finds 



v k +i -qk = v k - qk - 



(109) 



If this is combined with the induction hypothesis Vk-qk > 
qo • qk we obtain 



Vk+i qk>qo- qk- 



(110) 



The angle between qk+i and qo is larger than the angle 
between qk and qo, and these angles are in the range 
[0,7r]. According to Lemma ^ with q' = qk and q" = 
qk+i, it follows from Eq. I|110l) that 



Vk+i ■ qk+i > qo ■ qk+i- 



(ill) 



Since the induction hypothesis is true for k = 0, we ob- 
tain 



Vk ■ qk > qo ■ qk, 



(112) 



for all < k < M. This is true regardless of the choice 
of step size Ds. 

Suppose there exists a point s for which the inequality 
v(s) ■ q(s) > q(0) ■ q(s) does not hold. Then there exists 
a number a > 0, such that [v(s) — q{0)] ■ q(s) — — a < 0. 
By continuity of v and q, there must exist an interval 
si < s < S2 such that 

[v(s) - q(0)} ■ q(s) < -a/2, Vs G [si,s 2 ]. (113) 

Now consider approximations v^^s) in this interval. 

[ V W(s)-q(0)]-q(s) < [v™ (s) - v(s)} ■ q(s) - a/2 

< sup_ \\v^(s)-v(s)\\ 

se[si,s 2 ] 

-a/2. (114) 
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Since «( a )(s) converges uniformly to v(s), as 
Ds — ► 0, there exists a sufficiently small Ds such 
that sup||t> ( ' Ds) (s)-'y(s)|| < a/2. Thus, for a sufficiently 
small Ds it follows that [u( a )(s) - q(0)] • q(s) < 0, for 
all s 6 [Si,a^]< This contradicts Eq. I|112(l . since we may 
choose s to be some element of the partition (sk)^L 
inside the interval [si, S2]. (Such an element exists if Ds 
is sufficiently small.) Hence, by contradiction we have 
proved the statement of the lemma. □ 



1. Global search 

In the following we let Yp denote the solution of 
Eq. (|27[1 for a given run time T, with A = and B = 1. 
Similarly we let pp denote the solution of Eq. (|76|l . for 
the run time T . 

Lemma 3 

Y T {s)>Y {s), T>0, 0<s<l. (115) 

Proof. If p T (Q) = \E o (0))(E o (0)\, then, by construc- 
tion, we have 

Y T (s) = 2(Eo(s)\pt(s)\E (s)) - 1 = v T (s) ■ q(s), (116) 

where vt is the Bloch vector representing pp. Similarly 
Y (s) — Vq(s) ■ q(s), where v (s) is the Bloch vector 
corresponding to po(s) which is the solution of Eq. i|76|) 
with T = 0, and initial condition p (0) = \E o (0))(E o (0)\. 
The solution is po(s) = \E o (0))(E o (0)\. We may conclude 
that v (s) = q(0). Hence, Y T (s) > Y (s) is only Eq. 
in disguise. □ 

If we combine Eq. JSTJ) with Eq. (|115|l and let A = 0, 

B = 1, we find the following inequality 



1-F T (1)>J , 



(117) 



where 



e -nQ(s')-Q(s")] Z ( s >)Z(s")Y Q (s")ds"ds' 

(118) 

and where we have used that Yr(0) = 1. 

This implies that if lr(l) is to approach 1, then lo has 
to go to zero, or be negative, since Yr < 1. However, 
Iq > 0, as is shown below. 

Equation ||STJ yields 



Y o (s)=q(s)-q(0)> 



-(N-l)(2s-l) 



AT v /l + (AT-l)(2s-l) s 



(119) 



By inserting Eq. I|119|) into Eq. (|118|l and by making 
the change of variables x — \/N — l(2s' — 1) and y — 
y/N- l(2s" - 1) one obtains 



h > 



N 



N 



-ye 



-T[Q( K (x))-Q( K (y))} 



In the above expression we have introduced the function 
k(x) = 1/2 + x/(2y/N - 1). Furthermore, 

TQ(k(x)) = t£ (X) l-[l + (N-l)(2s'~l)rds' 
= a[$(z) - $(-\A/V - 1)], (121) 

where 



2N°y/N - 1 



and 



= / (1 + x ' 2 Ydx'. 



(122) 



(123) 



Note that $(-2;) = -$(x). By using Eq. 1121jl one can 
write 

TQ{k{x)) - TQ( K (y)) = a [$(1) - . (124) 
Define 

I(a,(3)= [ f a (x) [ g a {y)dydx, (125) 

J-/3 .7-/3 

From Eq. I|12U|) it follows that Iq fulfills the inequality 



Io>\l^-I(a,VN~l). 
The definition in Eq. I(125|) yields 



(126) 



—I(a,/3) = f a (/3) / g a (y)dy + g a (~(3) / f a (x)dx 
dp J-p J-p 

- Q *(/3) y S i n h(a$(y)) 

T+WJo (1 + y 2 ) 3 / 2 y 

^ e -"*(/3) /•£ cosh(a$(x)) , 
+2- „„,„,„ / = ax. 



(127) 



' ( 1 + ^2)3/2 J q l +x 2 



-vTv-T J-Vn-i 



(1+.T 2 )(1+ y 2 ) 3 /2 

(120) 



dydx. 



The second equality in the expression above is obtained 
by separating the integrals "J^V' into "J° „ + J Q ^" and 
make the change of variables x — ► —a; (and y — > — j/) in 
the J integrals. 

If one uses the inequality 

7; — ^v?^^^=t^t> Vye{O,0\. (128) 
(1 + y 2 ) 3 / 2 yiT^i + y 2 

in Eq. 1)127(1 ■ one obtains 

(129) 
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Since the integrand of the above integral is strictly posi- 
tive, one can conclude that 



— /(a,j8)>0, /3>0. 



(130) 



Hence, for a fixed a the function I (a, (3) is strictly in- 
creasing. Similarly, y/(N - 1)/TV increases with TV. If 
we assume TV > 2 it follows that 



Jo > ^I(a,VN=l)>-j=I(a,l) 
> ^=^F(a,0')df3'>O. 



(131) 



We can conclude that if Jo — > 0, then I(a, 1) — > 0, nec- 
essarily. This, in turn, implies F(a, /3')d{3' — > 0. We 
now investigate how the last expression depends on a. 

/- / F(a, 0)d(3> = f 1 ^F(a, (3')d(3\ (132) 
da J J da 



dF 
da 



(1+/? 2 ) 3 / 2 



-Q*(/3) 



(l + /?2)3/2 7 Q 



<5>(x)( 



-dx. 



(133) 



Both terms on the right-hand side in the above expression 
are strictly less than zero, for (3 > 0. It follows that 
Jp 1 F(a, &')d[3' is a strictly decreasing function of a. 

Thus, a necessary condition for the integral to go to 
zero is that a goes to infinity. Due to the form of a in 
Eq. (jnH it follows, for large TV, that TN-^ +1 / 2 1 -> oo 
is a necessary condition for ir(l) — *■ 1. Since ir(l) ~ * 1 
is equivalent to p — ► 1, we have shown that the suffi- 
cient condition T 3> A r<T+1 / 2 , proved in Sec. IV Al is also 
necessary. 



2. Local search 

The proof of the necessary condition in the local case 
resembles rather closely the proof of the global search 
case. However, we assume that a > 1. We let Yt denote 
the solution of Eq. Q27J) for a given run time T, with 
A = and B = 1, and with Q replaced by Q defined 
in Eq. (|STJ. Similarly we let pr denote the solution of 
Eq. H77|) . for the given run time T. 

It is straightforward to prove the counterpart of 
Lemma [3] 



Y T (r)>Y (r), T > 0, < r < 1. 



(134) 



Similarly we can obtain 1 — 3^(1) > Iq with Iq as in 
Eq. 1)1181) . but with Q replaced by Q. Moreover, 

-(N- l)(2r- 1) 



Y o (r) = q(r)-q(0)> 



N^l + (N - l)(2r- l) 2 



(135) 



Due to Eq. (|135H , Eq. (|120|) remains true if we replace Jo 
with Jq and Q(n(x)) with Q(k(x)). We also obtain 



TQ(k(x)) = 5[$(as) - 



where 



and 



T 



2LN°-WN- 1 2^^ arctan(VA r - 1) 



(l + a ; ' 2 ) CT - 1 ^ , . 



(136) 



(137) 



(138) 



Now we define I(a, /?) as in Eq. I|125|) but with a re- 
placed by S, and / Q , g a replaced by 



1 



1 



9a(y) 



— e a*W. 



y2)3/2 • 

(139) 

The line of reasoning from Eq. i|126|) to Eq. (| 1 33l> is unal- 
tered up to a replacement of s,a, Jo, and I(a,/3) with 
r, 5, $, J , and J(S, /?). We define F(a, (3) in analogy with 
F(a,f3) in Eq. 

Analogously to the global search case we find that a 
necessary condition for the limit Yt(1) — > 1 is that a — > 
oo, since J Q F(a,f3')dj3' is a strictly decreasing function 
in a. For large N, it follows, due to the form of a in 
Eq. (|137fl that TN~ a — > co is a necessary condition for 
p — > 1. Thus, we have proved that the sufficient condition 
for adiabaticity, proved in Sec. IV Bl is also necessary. 



VI. 



SUCCESS PROBABILITIES AND THE p 
LIMIT 



We have so far investigated the conditions for the suc- 
cess probability p to approach unity, which echoes the 
formulation of the standard adiabatic theorem. However, 
for an actual implementation of the adiabatic quantum 
computer we have to be satisfied with a success proba- 
bility less than unity in order to have a finite run time. 
It thus seems reasonable to ask how much of the results 
obtained for the p — > 1 limit survives when p < 1. 

Consider the set of pairs of list lengths and run times 
(TV, T) that gives precisely the success probability p. One 
may note that T does not necessarily have to be a func- 
tion of TV, as there may accidentally be several values 
of T, for a fixed TV, which give the same p. We show 
that the run time T is bounded by curves with the same 
asymptotic behavior as found in the p — > 1 cases. 

We assume that the system initially is in the ground 
state, i.e., poo(0) = 1 and poi(0) = 0. We denote 
P = Poo(l)- The strategy is to obtain expressions of 
the form 1 - p < U{T/N U ) or 1 - p > V{T/N"). If U 
and V are continuous and strictly decreasing functions, 
we may invert them, on suitable domains, and obtain 
bounds N v V- x (l-p) < T < N v U~ x (l-p). We derive 
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both upper and lower bounds in the wide-open case. In 
the semi-open case we obtain upper bounds only. 

One may note that the upper bounds can be seen as 
sufficient conditions on the run time to achieve a success 
probability which is at least p. 



A. The wide-open case, global and local 

An upper bound for the global wide-open case can be 
obtained directly from Eq. I|8()|l. The result is 



T < N cr+1 / 2 - 



1-p' 



a > 0. 



(140) 



We obtain a lower bound by using results from Sec. IV C 11 
By combining Eqs. I|117|) and i|131|) one obtains 



1-p > 



2V2Jo 



F(a,/3)d/3 = C(a), 



(141) 



where F is defi ned in Eq. !|129|) . and a in Eq. (|122|l . Wc 
know from Sec. IV C"T1 that C is strictly decreasing in a. 
Moreover, C is continuous. Hence, C is invertible and 



T > 2AVA - 1C _1 (1 - p), N>2. 



(142) 



The bounds in Eqs. 1140(1 and (|142l) fulfill the same 
jyer+1/2 ra ^ e ag we f ounc [ j n the p — ► 1 investigation. 

For the wide-open local case we may use Eq. I)84|l . and 

P = Poo(l) = poo(l) to obtain 



T < 



TT 2 N a 

Yi -p 



(143) 



Analogously to the global search case we can obtain 
a lower bound to the wide-open local search from 
Sec. EIH as 



1-P > 



1 



2V2 JO 



F(a 1 P)dp = C(a), 



(144) 



with a as in Eq. (|137[) . This results in the bound 

T > 2A CT arctan( v / A T7 T)C*- 1 (l -p), (145) 

which, for large N, increases as N" . 

We have found that for a fixed success probability p the 
run time is locked between two bounds. These bounds 
have the same scaling with the list length N as found in 
the p — > 1 limit. 

B. The semi-open case, global and local 

In the semi-open global search case we use Eq. ()54|l . 
with poi (0) = and poo(0) = 1 to obtain 



T < -kN 



BK + 6A 
A 2 (l-p)' 



(146) 



In the semi-open local search case one may use Eq. <|72|) 
to obtain 

The bounds in Eqs. (|14(j[) and 1)1471) are sufficient con- 
ditions on the run time to obtain at least the success 
probability p, for a given list length N ■ 



C. Numerical analysis 




log 2 N 



FIG. 1: (Color online) Global search with W(s) = H(s). 
Each curve shows log 2 T vs log 2 N, where T is the run time 
needed to obtain a given success probability p at s = 1, and 
where iV is the list length. The dimensionless run time T is 
the quotient between the physical run time and To defined in 
Eq. The dotted curves correspond to the wide-open 

case, defined by the degree of openness u> — 1. The dashed 
curves correspond to the semi-open case ui = 0.9, and the 
solid correspond to u> = 0.5. Within each group, each curve 
corresponds to a success probability, which counted from be- 
low is p — 0.4, 0.5, 0.6, 0.7, 0.8. As seen, the asymptotic slope 
appears to be independent of the choice of success probability 
p. For the wide-open case the asymptotic slope seems to tend 
to v — 3/2, independent of p. For the semi-open cases the 
asymptotic slope appears to tend to v — 1, independent of p. 

We have numerically calculated the actual run times 
of the wide-open global and local search. Let us choose 
A = cos(u;7r/2) and B = sin(ciJ7r/2). This choice of A and 
B is merely a convenient way to obtain an interpolation 
between the closed case u> — and the wide-open case 
w = 1, without effectively changing T. Moreover we have 
chosen W(s) — H(s). Thus, we consider an environment 
that monitors the energy of the system. The initial con- 
dition at s = is p(0) = \ip)(ip\ with |i/>) as in Eq. JUJl. 
In other words, we begin at the initial ground state of 
the Hamiltonian in Eq. (|13f) . 

We should also comment upon the units. As seen in 
Sec. IIIII we have put the energy gap between the ground 
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state and the degenerate first excited state to 1 in the 
initial Hamiltonian H(0). Depending on the actual sys- 
tem used to realize the search procedure, the energy gap 
is actually AS in some appropriate energy unit. To put 
the energy gap to 1 corresponds to measuring energy in 
units of A£. Furthermore, as we have chosen h = 1, it 
follows that the run time T, used throughout this paper, 
can be seen as the quotient T = T/Tq, of the physical 
run time T and To defined by 

T = A£ _1 . (148) 

In Fig. ^ the dotted curves correspond to the wide- 
open global search with success probability ranging from 
p = 0.4 to p = 0.8. The curves seem to follow closely to 
a linear asymptote with slope v = 3/2, independent of 
the choice of p < 1. The dotted curves in Fig. |5] similarly 
correspond the wide-open local search case. The asymp- 
totic slope of these curves seems to be v = 1 for every 
choice of p < 1 . 

We have also calculated the actual run times of the 
semi-open global and local search numerically. Again we 
use W(s) — H(s), A = cos(w7r/2), and B — sin(w7r/2). 
The results for the global case are shown in Fig.(2 in form 
of the dashed curves corresponding to u> = 0.9, and the 
solid curves corresponding to lu — 0.5. All these curves 
seem to tend to the asymptotic slope v = 1, which is in 
line with our analytical results. In Fig. [21 we similarly 
have created plots for the semi-open local search case. 
These again correspond to the dashed curves with u) = 
0.9, and the solid curves with uj = 0.5. The results are 
qualitatively similar to the semi-open global search case, 
but with asymptotic slope v = 1/2. 

Figures |21 and 0] show the global and local search case 
for a fixed success probability p = 1/2 but for different 
degrees of openness uj. Figure|3suggests that in the semi- 
open cases the asymptotic slope is v = 1, while in the 
wide-open case the asymptotic slope is v — 3/2, which is 
consistent with our analytical results. In Fig.^Jthe curves 
in the local search case appears to have a qualitatively 
similar behavior to the global search case, except that 
the asymptotic slope is v — 1 in the wide-open case, and 
v = 1/2 in the semi-open case. 

From Figs. |3] and 0] one may obtain an intuitive un- 
derstanding of the "discontinuity" of the scaling of the 
run times, when going from the semi-open to the wide- 
open case. The closer we get to the wide-open case (i.e., 
the closer uj gets to 1), the higher up in list-lengths the 
curve "follows" the curve of the wide-open case, before it 
turns to its asymptotic slope. In this context one should 
also note that for every fixed list-length N, the system is 
continuous with respect to u>. In other words, the discon- 
tinuity only appears for the asymptotes in the N — ► oo 
limit, not in any system corresponding to a fixed list. 
The list-length N corresponds to the dimension of the 
total Hilbert space, but in the effective two-dimensional 
Hilbert space the list-length only enters as a parameter 
in the Hamiltonian. If we are to understand the dis- 
continuity of the asymptotic scaling as a function of the 
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FIG. 2: (Color online) Local search with W(s) — H(s). Each 
curve shows log 2 T vs log 2 N, where T is the run time needed 
to obtain a given success probability p at s = 1, and where 
N is the list length. The dimensionless run time T is the 
quotient between the physical run time and To defined in 
Eq. 11481 . Similar to Fig. Q the dotted curves correspond 
to the wide-open case, defined by the degree of openness 
ill = 1. The dashed curves correspond to the semi-open case 
uj = 0.9, and the solid curves correspond to ui — 0.5. Within 
each group each curve corresponds to a success probability 
p = 0.4, 0.5, 0.6, 0.7, 0.8, counted from below. Similarly to the 
global search, the asymptotic slope appears to be independent 
of the choice of success probability p. For the wide-open case 
the asymptotic slope seems to tend to v = 1, independent of 
p. For the semi-open cases the asymptotic slope appears to 
tend to the v — 1/2, independent of p. 



degree of openness, in the case W(s) — H{s), we may 
consider how the parameter N enters the terms A[H,p] 
and B[H, [H, p]\ in the master equation. The minimum 
energy gap of the Hamiltonian is proportional to l/-y/iV. 
Thus, at the minimum gap the list length TV enters as 
A/^/N in the Hamiltonian part of the master equation, 
but as B/N in the decohering part. No matter how small 
A is, as long as it is nonzero, A/y/N always dominates 
B/N in the limit of large N. Although this, strictly 
speaking, only holds at s = 1/2, it seems intuitively rea- 
sonable that the Hamiltonian term dominates the evolu- 
tion, for large N, as long as A ^ 0. Intuitively it also 
seems reasonable that if the Hamiltonian dominates the 
evolution, then the Hamiltonian determines the asymp- 
totic scaling of the run time, and thus the scaling becomes 
the one of the ideal closed case. However, if A — then 
the decohering term determines the scaling. 



VII. CONCLUSIONS 

We investigate the effect of decoherence on the adi- 
abatic search of a marked element in a disordered list 
of N elements. More specifically we consider deco- 
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FIG. 3: (Color online) Global search with W(s) = H(s) 
and success probability p — 1/2. The curves show log 2 T vs 
log 2 N, where T is the run time needed to obtain the suc- 
cess probability p = l/2ats=l, and where N is the list 
length. The dimensionless run time T is the quotient be- 
tween the physical run time and To defined in Eq. (11481 . 
Each line shows the result for a given value of the degree 
of openness uj. Counted from below, the lines correspond to 
lj = 0, 0.1, . . . , 0.9, 1. Hence, the uppermost line corresponds 
to the wide-open case u) = 1 and the lowermost line corre- 
sponds to the closed case ui — 0. All the lines seem to tend to 
the slope v = 1, except the uppermost which seems to tend 
to v = 3/2. 



herence with respect to the instantaneous eigenbasis of 
the adiabatic search Hamiltonian. The decoherence is 
modeled by using a master equation of the form p = 
-iA[H,p] - B[W, [W,p\] (with h = 1), where H is the 
time-dependent search Hamiltonian and W is a time- 
dependent Hermitian operator that commutes with H. 
The operator W is only subjected to the condition in 
Eq. which, in some sense, states that the fluctua- 

tions of W (with respect to time) should not grow with 
the list-length N. 

We prove sufficient conditions for the p — ► 1 limit of 
the success probability p. In the semi-open case (A ^ 0) 
we prove that these sufficient conditions are T ^> N and 
T » >/N in the global and local search case, respectively, 
where T is the run time of the adiabatic search. In other 
words, the results for the ideal adiabatic search (B = 0) 
remains valid even in the presence of decoherence with 
respect to the instantaneous eigenbasis. 

In the wide-open case (A = 0), sufficient and necessary 
conditions for the p — > 1 limit, are deduced. In the case 
where W — H we find that the conditions for p — > 1 are 
T > 7V 3/2 and T > N for the global and local search, 
respectively. 

We finally investigate how the run time T depends on 
N for a fixed success probability p less than unity. We 
find that it is sufficient with run times that scales like N 
and >/N, in the semi-open case, as for the ideal global 
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FIG. 4: (Color online) Local search with W(r) = H(r) 
and success probability p = 1/2. The curves show log 2 T 
vs log 2 N plots, where T is the run time needed to obtain 
the success probability p = 1/2 at s = 1, and where N is 
the list length. The dimensionless run time T is the quotient 
between the physical run time and To defined in Eq. 11481 . 
Each line shows the result for a given value of the degree of 
openness uo. Counted from below, the lines correspond to 
uj = 0, 0.1, . . . , 0.9, 1. All the lines seem to tend to the slope 
v = 1/2, except the uppermost which appears to tend to 
v = 1. 



and local adiabatic searches. In the wide-open case we 
instead find the scalings N 3 / 2 and N, for the global and 
local search, respectively. 

In view of these results the following picture emerges. 
In the semi-open cases the Hamiltonian dynamics dom- 
inates the behavior of the system for sufficiently large 
list lengths. To be more precise, an increased degree 
of eigenbasis decoherence does increase the run time of 
the adiabatic search, but only through the constant C 
of the asymptote CN V . In the wide-open case, however, 
the protective effect of the Hamiltonian dynamics is not 
present. Consequently, the asymptotic behavior depends 
directly on the choice of W. As a consequence the asymp- 
totic behavior of the wide-open case may be distinctly dif- 
ferent from that of the semi-open case, although we have 
the same choice of W. One may note that the abrupt 
change in asymptotic behavior with respect to the tran- 
sition from semi-open to wide-open system, raises a warn- 
ing sign concerning the use of the wide-open evolution as 
an approximation to an "almost" wide-open evolution. It 
should be noted, however, that for any fixed list length, 
the system is continuous with respect to the degree of 
openness, and that the discontinuity rather concerns the 
asymptotic scaling of the run time in the limit of large 
list-lengths. 

The decoherence model provided by Eq. implies 
that we disregard the memory effects in the environment 
and the corresponding fluctuations, which seems justifi- 
able if the relevant time scales of the evolution is much 
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larger than the dominating time scales of the memory of 
the environment. The modeling of memory effects gen- 
erally leads to integro-differential equations [53, 13 i or 
to time-local non-Markovian equations [2{|. To include 
such effects into the analysis of the influence of the envi- 
ronment on the efficiency of the adiabatic quantum com- 
puter would be interesting but challenging. 

Finally we note that, in spite of the apparent differ- 
ence between the adiabatic quantum computing scheme 



and the traditional circuit model, it has been shown 
that these two models are, in a certain sense, equivalent 
[lEIlSl- However, this equivalence does not concern the 
robustness to noise, relaxation, or decoherence. Previous 
findings suggest that the adiabatic quantum computer to 
a certain extent should be resilient against various kinds 
of open system effects. The results of this paper provide 
further evidence for this. 
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